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ABSTRACT 

We present the results of a study of a statistically significant sample of galaxies which 
clearly demonstrate that supermassive black holes are generically present in all mor- 
phological types. Our analysis is based on the quantitative morphological classification 
of 1.12 million galaxies in the SDSS DR7 and on the detection of black hole activity 
via two different methods, the first one based on their X-ray/radio emission and the 
second one based on their mid-infrared colors. The results of the first analysis confirm 
the correlation between black hole and total stellar mass for 8 galaxies and includes 
one galaxy classified as bulgeless. The results of our second analysis, consisting of 
15,991 galaxies, show that galaxies hosting a supermassive black hole follow the same 
morphological distribution as the general population of galaxies in the same redshift 
range. In particular, the fraction of bulgeless galaxies, 1,450 galaxies or 9%, is found 
to be the same as in the general population. We also present the correlation between 
black hole and total stellar mass for 6,247 of these galaxies. Importantly, whereas pre- 
vious studies were limited to primarily bulge-dominated systems, our study confirms 
this relationship to all morphological types, in particular, to 530 bulgeless galaxies. 
Our results indicate that the true correlation that exists for supermassive black holes 
and their host galaxies is between the black hole mass and the total stellar mass of 
the galaxy and hence, we conclude that the previous assumption that the black hole 
mass is correlated with the bulge mass is only approximately correct. 

Key words: galaxies: general — galaxies: Seyfert — galaxies: active — galaxies: 
spiral — galaxies: bulges — galaxies: evolution — infrared: galaxies 
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1 INTRODUCTION 

According to the standard picture, the structure we observe 
today in our Universe formed as a consequence of the growth 
of fluctuations in the primordial dark matter distribution 
which gravitationally attracted gaseous baryonic matter 
that later settled into disk galaxies (e.g. White & Rees 
1978; White & Frenk 1991). Galaxy mergers within halos 
subsequently transformed disks into bulges and the gas fell 
to the center, triggering starbursts and feeding the rapid 
growth of black holes (BH). These in turn responded by feed- 
ing energy back to the surrounding gas (e.g. Sanders et al. 
1988; Hopkins et al. 2008). This picture of structure for- 
mation is based on observational work which, over the past 
few decades, indicates that most -if not all- massive galax- 
ies with a spheroidal component have a supermassive black 
hole (SMBH) at their center (e.g. Kormendy & Richstone 
1995; Magorrian et al. 1998). In addition, estimates of the 
BH masses of these giant galaxies have been found to tightly 
correlate with the spheroid/host luminosity (or stellar mass) 
and also with the stellar velocity dispersion within the 



host bulges (e.g. Gebhardt et al. 2000; Ferrarese & Merritt 
2000). 

However, there is growing evidence that SMBH not only 
exist at the centers of massive galaxies with a central bulge, 
but also at the centers of galaxies with pseudo-bulges and 
that these systems also exhibit a correlation between black 
hole mass and total stellar mass (Simmons et al. 2011). 
This supports the alternative view that SMBH growth 
might be due to minor mergers and/or secular processes 
(e.g. Green et al. 2010; Simmons et al. 2011). Although, 
as pointed out by Okamoto (2012), the main channel of 
pseudo-bulge formation may in fact not be secular disk evo- 
lution but rather rapid gas supply at high-redshift. 

Given the uncertainties concerning the current under- 
standing of the formation of bulges and pseudo-bulges in 
galaxies, the problem of identifying the main process at 
play in the formation of SMBH at the centers of galaxies 
with a spheroidal component remains a difficult challenge at 
present. On the other hand, the recent discovery of a SMBH 
at the center of Henize 2-10 (Reines et al. 2011), a dwarf 



© 0000 RAS 



2 Marleau et al. 



starburst galaxy lacking any substantial spheroidal compo- 
nent, suggests that a merger is not essential to making a 
SMBH and that BH growth may after all require substan- 
tial and ongoing secular processes, at least in these galaxies. 
Moreover, if SMBH are also generic to bulgeless galaxies 
and dwarf galaxies, as indeed the current work will present 
evidence to demonstrate they are, these galaxy types may 
offer the best probe of SMBH formation and evolution in 
the local Universe. In this regard, a SMBH hosted by a 
bulgeless galaxy and/or dwarf galaxy is more likely "pris- 
tine", as it is unlikely that it has experienced mergers and 
dynamical interactions but has rather undergone a quieter 
evolution. Moreover, pure-disk galaxies are far from rare: lo- 
cally, only four of the nineteen giant galaxies are ellipticals or 
have classical bulges (Kormendy et al. 2010). While in the 
Sloan Digital Sky Survey (SDSS), based on the quantitative 
morphological analysis of 1.12 million galaxies in the DR7 
(Simard et al. 2011, hereafter Sll), the fraction of galax- 
ies with an optical Sersic profile index rig < 1.5 and bulge 
fraction B /T < 0.1, i.e. nearly bulgeless, is ~14%. How- 
ever, despite the large fraction of bulgeless galaxies and their 
clear importance in challenging our current understanding 
of galaxy formation and evolution, a systematic survey of 
SMBHs in bulgeless galaxies has never been undertaken, 
though this is in part due to the lack of a well-defined sam- 
ple. In this paper, we attempt to rectify this situation and 
substantially extend the status of the present knowledge of 
SMBH and their host galaxies by presenting the first sys- 
tematic survey of SMBHs in galaxies of all morphological 
types, i.e. including those classified as bulgeless galaxies, 
based on the quantitative morphological analysis of 1.12 mil- 
lion galaxies in the SDSS DR7 (Sll). 

In this work we report the results of two separate studies 
which explore the nature of SMBH and their host galaxies, 
the first based on an X-ray/radio selected sample and the 
second on an infrared selected sample. The structure of our 
paper is as follows. In Section 2, we describe the method that 
we used to select our sample of morphologically classified 
galaxies. In Section 3, we present the selection and analysis 
of our X-ray/radio selected sample of active SMBH, while 
in Section 4, we present the selection and analysis of our in- 
frared selected sample of active SMBH. In Section 5, we sum- 
marize our results. Considered together, we believe these two 
studies are complementary in nature, in that while our first 
analysis of an X-ray selected sample provides a well-tested 
and robust estimate of black hole mass but with a small 
statistic, our second analysis of an infrared selected sample 
provides estimates on the black hole mass, parametrised by 
their accretion rate, but with a large statistic. 
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Figure 1. Rcdshift distribution of the 688 galaxies in the CSC- 
SDSS catalog with structural parameters. 



galaxies. In addition to these two criteria, the sum of the 
flags DEBLENDED_AS_PSF and SATURATED were re- 
quired to be zero to eliminate objects that were found to 
be unresolved children of their parents as well as saturated 
objects. The total number of objects in the Sll catalog is 
1,123,718. 

Galaxy structural parameters obtained by Sll were 
measured from bulge -|- disk decompositions performed 
using version 3.2 of the GIM2D software package 
(Simard et al. 2002). In addition to the canonical de Vau- 
couleurs bulge (Sersic index nt = 4) -I- exponential disk 
fitting model (Sll, Table 1), two other fitting models were 
used. The first additional model was a free nj, bulge + disk 
model, with nt allowed to vary from 0.5 to 8 (Sll, Ta- 
ble 2), and the second one was a single component pure 
Sersic model, with Ug also allowed to vary from 0.5 to 8 
(Sll, Table 3). 

In order to classify a galaxy as bulgeless, we required 
an optical Sersic profile index Ug < 1.5 and a bulge fraction 
B/T < 0.1. For the whole Sll catalog, this yielded a list 
of 88,594 galaxies, corresponding to about 14% of the total 
sample. 



2 QUANTITATIVE MORPHOLOGICAL 
PARAMETERS 

As the foundation of all our work we used the quantitative 
morphological classification of Sll to extract structural pa- 
rameters for 1.12 million galaxies. In this work, the sources 
that were fitted were selected from the Legacy area with 
14.0 ^ rnpetro,r,corr ^ 18.0, where mpetro,r,corr IS the r-baud 
Petrosian magnitude corrected for Galactic extinction ac- 
cording to the extinction values given in the SDSS database, 
and morphological type Type — 3, i.e. sources classified as 



3 X-RAY SMBH SAMPLE SELECTION 

3.1 CSC-SDSS Cross-Match Catalog 

In our first analysis, which we now describe, we investigated 
the properties of a sample of galaxies containing AGN se- 
lected via their X-ray detection. In order to achieve this 
purpose we utilized the Chandra Source Catalog/Sloan Dig- 
ital Sky Survey (CSC-SDSS) cross-match catalog to select 
active black holes via their X-ray emission. The X-ray color 
(or hardness ratio), often used to identify AGN, was left un- 
constrained in this selection process in order to avoid filter- 
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Figure 2. Distribution of the Sersic index rig (left), the bulge fraction {B/T)r {middle) and the bulge Sersic index nj, (right) of the 
galaxies in the CSC-SDSS catalog with structural parameters. 



ing out some AGN candidates, as, for example, many of the 
brightest X-ray sources are identified as AGN candidates by 
WISE, regardless of their hardness ratio (Wang et al. 2004; 
Stern et al. 2012). In total, 8,997 galaxies in the SDSS were 
cross-matched to a Chandra X-ray detection. Of these, we 
subsequently identified 688 as having been assigned struc- 
tural parameters in the catalog of Sll. Finally, after apply- 
ing our selection criteria of rig < 1.5 and B/T < 0.1, we 
were left with a sample of 26 galaxies classed as bulgeless. 



3.2 Redshift Distribution and Structural 
Parameters 

In order to gain some insight into the nature of our of X-ray 
selected sample of 688 galaxies with structural parameters, 
we examined their redshift and structural parameter distri- 
butions and compared them to the general population. The 
redshift distribution of the X-ray selected sample has been 
plotted in Figure 1, while the distributions in Sersic index, 
bulge fraction and bulge Sersic index are depicted in Fig- 
ure 2. The Sersic profile index was plotted only for galaxies 
better fitted with a Sersic model, i.e. with a F-test probabil- 
ity Pps > 0.32 (Sll), whereas both the bulge fraction and 
the bulge Sersic index were plotted only for galaxies better 
fitted with a bulge + disk components model, in this case 
with Pps ^ 0.32. The shape of the redshift distribution was 
found to be similar to that of Sll (see their Figure 2), with 
a peak at z ~ 0.1, indicating no obvious redshift bias in 
the X-ray sample. However, the distribution in Sersic profile 
was found to differ, being flatter at low rig, and possessing 
no prominent preferred value at Tig — 1, contrary to what 
was seen in the general population (see Figure 6). Similarly, 
the distribution in bulge fraction was also found to show a 
deficit of {B/T)r < 0.3 as compared to the general popula- 
tion (see Figure 6). Taken together, therefore, the structural 
parameter distributions seemed to indicate a bias towards 
non-bulgeless galaxies in the X-ray selected sample. 



3.3 Stellar Mass Estimates 

One of the primary objectives of our work was to investigate 
the relationship between the stellar masses of our galaxies 
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Figure 3. Total stellar mass versus black hole mass estimated 
from the fundamental plane of black hole activity for our X- 
ray/radio sample (black filled circles). The bulgeless galaxy is 
shown as a red filled circle. For comparison, we also plot the em- 
pirical relation (dashed line) and the observed scatter of 0.3 dex 
(grey shaded region) of Haring & Rix (2004). 



and their supermassive black holes. In order to obtain esti- 
mates of the stellar masses of our sample, we made use of 
the catalog of Mendel et al. (2012, hereafter M12). In this 
work, M12 estimated the stellar masses of 669,634 galaxies 
in the Legacy area of the SDSS DR7, with spectroscopic red- 
shifts in the range 0.005 z 0.4. Importantly for us, M12 
considered only sources classified as galaxies both photomet- 
rically and spectroscopically by the SDSS, i.e. sources with 
photoObj.type = 3 and specObj.specClass = 2. A key conse- 
quence of the nature of this selection was that it should have 
had the effect of filtering out most of the unobscured (type 
1) AGN, i.e. those showing broad emission lines. Crucially, 
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this implies that our final sample possessing stellar masses, 
obtained from the catalog of M12, both for this analysis and 
the analysis in Section 4, should have been composed mainly 
of obscured (type 2) AGN. Although, this might be viewed 
as a restriction, as ideally one would like to have a sample 
composed of both types of AGN, there are two immediate 
benefits to having a sample composed of type 2 AGN. The 
first is that this provides an opportunity to study SMBHs in 
type 2 AGN on a large scale. The second benefit of obtaining 
such a sample is that one can expect a minimal contribution 
from obscured AGN to the total optical light in these galax- 
ies, and hence, a negligible effect from the AGN optical light 
to the estimates of their stellar masses. 

The stellar masses obtained by Ml 2 were derived via 
the fit of the spectral energy distribution (SED) of SDSS 
ugri photometry with the flexible stellar population synthe- 
sis code of Conroy et al. (2009). For our selected sample of 
bulgeless galaxies, we used the total stellar mass estimates 
from the Sersic profile fitting (M12, Table 2). While for the 
remainder of our sample, we used the total, bulge and disk 
stellar mass estimates from the de Vaucouleurs bulge -|- ex- 
ponential disk fitting (M12, Table 3). It is important to note 
that the fit to the total, bulge and disk photometry was car- 
ried out independently and, therefore, it is not necessarily 
true that M*;, + M*d = M^t. Of the 688 galaxies in our X- 
ray /morphology sample, we were able to match 404 galaxies 
with the M12 catalog and get their stellar masses. 
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Figure 4. Rcdshift distribution of the 15,991 WISE color selected 
AGN with structural parameters. 



Table 1. Stellar Mass Estimates 



3.4 Black Hole Mass Estimates from Fundamental 
Plane of BH Activity 

To compute an estimate of the black hole mass for our X-ray 
sample, we appealed to the fundamental plane of black hole 
activity (Merloni et al. 2003), which correlates the X-ray 
and radio luminosities with the mass of the black hole. 

In order to search for the level of radio emission from the 
detected Chandra X-ray sources in our galaxy sample, we 
utilized the source catalog produced by the National Radio 
Observatory/ Very Large Array (NRAO/VLA) Sky Survey 
(NVSS). This radio survey was undertaken by the NRAO 
VLA telescope and covered the sky north of a declination of - 
40 degrees, at a frequency of 1.4 GHz, a resolution of 45" and 
a limiting peak source brightness of about 2.5 mjy/beam. Of 
the 404 sources that we had previously identified as having 
been assigned a stellar mass in the M12 catalog, we man- 
aged to find 8 that also had an NVSS radio detection, one 
of which we identified as bulgeless. The stellar masses of 
our final 8 galaxies are given in Table 1, while their mor- 
phological parameters are listed in Table 2. We note that 
the galaxy 587724199889666262 can also be found in the 
2MASS-selected Flat Galaxy Catalog of Mitronova et al. 
(2004). 

The empirical relationship relating black hole mass to 
the emitted compact radio and hard X-ray luminosities, 
spanning nine orders of magnitude in black hole mass, is 
formalized by the equation 



log Lr = 0.60 logLx + 0.78 logM 7.33, 



(1) 



where Lh is the radio luminosity at 5 GHz in erg/s, Lx is 
the 2-10 keV X-ray luminosity in erg/s, and M is the mass 
of the black hole in solar masses. 
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Figure 5. Distribution of the Sersic profile index Ug (left), the bulge fraction {B/T)r {middle) and the bulge Sersic index (right) of 
the WISE color selected AGN with structural parameters. 
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Figure 6. Distribution of the Sersic profile index rig (left), the bulge fraction {B/T)r {middle) and the bulge Sersic index nf, (right) of 
the SDSS DR7 galaxies in the Sll morphological catalog. 



Table 3. Estimated Black Hole Mass 
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40.916 


38.634 


8.6598 


588015508739195010 


40.476 


38.924 


9.3690 


587731513142673426 


41.225 


38.999 


8.8898 


587731513142673623 


39.609 


37.689 


8.4525 


587731513678561426 


42.095 


38.634 


7.7521 


587731513687670913 


42.097 


39.146 


8.4072 


587724199889666262 


40.506 


37.871 


7.9969 


587730775500128413 


43.188 


39.427 


7.9291 



lar mass for our sample of 8 galaxies, including the one 
classified as bulgeless, is shown in Figure 3. For compari- 
son, we show in the same figure the empirical relation of 
Haring & Rix (2004) (dashed line) derived using dynami- 
cal and stellar masses of bulge- dommated systems, for which 
M*j, ~ Mtt. In contrast, our total stellar masses are based 
either on the Sersic profile fitting (for our bulgeless galaxy) 
or the de Vaucouleurs bulge -I- exponential disk profile fit- 
ting (for the remainder of the sample). The fact that the 
correlation between black hole mass and total stellar mass 
is consistent with the relation of Haring & Rix (2004) sug- 



gests that the true correlation might not be with the bulge 
mass but rather with the total stellar mass. That this is in- 
deed the case becomes clearer in our second study, which we 
now detail. 

4 INFRARED SMBH SAMPLE SELECTION 

4.1 WISE-SDSS Cross-Match Catalog 

In our second analysis, we utilized the Wide-field Infrared 
Survey Explorer (WISE) all-sky catalog to identify both un- 
obscured (type 1) and obscured (type 2) AGN. In this re- 
gard, it is important to note that most surveys of AGN that 
have previously been undertaken have been biased towards 
observing mainly unobscured (type 1) AGN, even though 
theoretical models generally predict a large population of 
obscured (type 2) AGN, outnumbering type 1 AGN by a 
factor of ~3 (e.g. Comastri et al. 1995; Treister et al. 2004; 
Ballantyne et al. 2011). A key advantage to using the mid- 
infrared selection of AGN is that it does not sulfer from 
this bias, as it relies upon distinguishing the approximately 
power-law AGN spectrum from the black body stellar spec- 
trum of galaxies which peaks at rest-frame 1.6^m, and is 
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Figure 7. Distribution of thic total stellar masses {left) and lower limits on black hole masses (Lboi = L^dd) fight) of the 6,247 WISE 
eolor seleeted AGN with structural parameters and stellar masses. 



thereby able to identify both unobscured (type 1) and ob- 
scured (type 2) AGN. Another advantage to working with 
mid-infrared data is that they allow AGN to be easily dis- 
tinguished from stars and galaxies. Furthermore, an added 
benefit of mid-infrared selection is that it does not suffer 
from dust extinction and is sensitive to the highest redshift 
sources. 

In order to obtain our infrared sample we downloaded 
part22 to partSO of the WISE All-Sky Catalog (covering the 
area covered by SDSS, i.e. dec > -23.1899 degrees). The 
WISE All-Sky Release Source Catalog contains positions 
and photometry at 3.4, 4.6, 12 and 22 micron for 563,921,584 
point-like and resolved objects detected on the Atlas Inten- 
sity images. Photometry was performed using point source 
profile-fitting and multi-aperture photometry and the esti- 
mated sensitivities are 0.068, 0.098, 0.86 and 5.4 mjy (5(t) 
at 3.4, 4.6, 12 and 22 micron in unconfused regions on the 
ecliptic plane. J2000 positions and uncertainties were recon- 
structed using the 2MASS Point Source Catalog as astro- 
metric reference. Astrometric accuracy is approximately 0.2 
arcsec root-mean-square on each axis with respect to the 
2MASS reference frame for sources with signal-to-noise ra- 
tio greater than forty. 

Based on their distinctive spectral energy distribu- 
tion in the infrared (Wu et al. 2012; Stern et al. 2012; 
Asset et al. 2010), we selected AGN in the WISE all-sky 
catalog out to a redshift of 2 ~ 3.5 by applying the mid- 
infrared color criterion Wl-W2^ 0.57 (Stern et al. 2012; 
Wu et al. 2012). We found that ~15% of the WISE sources 
satisfied this criterion (from a total of 3.27072 x 10* sources, 
we were left with 51,048,962.) 

Having obtained our initial sample we then used 
the US Virtual Astronomical Observatory (VAO) Cross- 
Comparison Tool to upload the WISE-selected AGN catalog 
and cross-matched this with the online SDSS DR 7 catalog, 



using a 0.5 arcsec match radius (Stern et al. 2012, posi- 
tional accuracy with 2MASS). The number of cross-matched 
sources equaled 6,969,533, corresponding to about 2% of the 
initial sample. In total, 15,991 galaxies were found to have 
structural parameters derived by Sll and we subsequently 
found that 1,450 of these galaxies (or about 9%) satisfied 
our bulgeless selection criteria. 

4.2 Redshift Distribution and Structural 
Parameters 

The distributions of our WISE color selected AGN galaxy 
population have been plotted in Figures 4 and 5, as a func- 
tion of their redshift, z, and structural parameters, Ug, 
{B/T)r and ni,, respectively. As in Figure 2 of Section 3.2, 
the Sersic profile index was plotted only for galaxies with 
PpS > 0.32 (Sll), whereas both the bulge fraction and 
the bulge Sersic index were plotted only for galaxies with 
PpS ^ 0.32. Simple consideration of these distributions lead 
to some important new results regarding the properties of 
AGN and, therefore, more generally for galaxies containing 
SMBH. 

The first, and most important of these statistical re- 
sults, which can be immediately seen in the distribution of 
the Sersic index of our infrared selected sample. Figure 5 
(left), is that the AGN are smoothly and significantly dis- 
tributed over all values of Tig. This necessarily implies that 
AGN in our sample, and therefore their associated SMBH, 
are present in galaxies of all morphological types. This is a 
key discovery, as until now relatively few SMBH were known 
to exist in galaxies that did not possess a significant bulge 
component. 

The second important statistical result, is that the mor- 
phological distributions of the AGN sample and the general 
population of galaxies have a very similar form. This can be 
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Figure 8. Left: Total stellar mass versus black hole mass obtained using the bolometric luminosity for our infrared sample of 5,718 
galaxies not classified as bulgeless {black open circles) and for the 545 bulgeless galaxies {green filled circles). The data points and 
bisector linear regression fit to the data {dashed line) are plotted for Lhoi/LEdd = 0.1. Also shown are the fit for Li,oi/ L^dd = 1-0 {dotted 
line) and for Lf^^i/ L^^^ = 0.01 {solid line). For comparison, we also plot the observed scatter of 0.3 dcx of the empirical relation of 
Haring & Rix (2004) {grey shaded region). Right: Same as left figure, showing the galaxies with 3.9 < nj, < 4.1 highlighted as red filled 
circles. 



clearly seen by comparing the Sersic distribution of the in- 
frared selected AGN, Figure 5 (left), with that of the parent 
population Figure 6 (left). Alternatively, comparing Figure 5 
(middle and right) with Figure 6 (middle and right), simi- 
larly reveals that in the case when the galaxies have been 
fitted with a bulge + disk component model, the galaxies 
with a SMBH again follow a very similar distribution of 
bulge fraction and bulge Sersic index to the general popu- 
lation of galaxies. Hence, from these comparisons we learn 
that not only are the SMBH in our infrared selected sam- 
ple distributed over all morphological galaxy types, but that 
the form of their morphological distribution is very similar 
to that of the general population of galaxies, over the same 
redshift range, i.e. 0.05 < z < 0.25. Moreover, this similar- 
ity property is essentially robust with respect to the form 
of morphological fitting that is used. An immediate conse- 
quence of this approximate scaling relation between the dis- 
tributions of the AGN and the general population, is that it 
implies that the fraction of galaxies that contain a SMBH is 
approximately the same for each morphological type. Alter- 
natively, one may say that the fraction of a particular type 
of galaxy is the same in both the AGN and general popu- 
lations. In particular, the fraction of bulgeless galaxies that 
host a SMBH is the same as the fraction of bulgeless galax- 
ies in the general population, which is ~9% in the redshift 
range considered. 

An important set of questions concern the extent to 
which our statistical sample may be assumed to be repre- 
sentative of the entire population of AGN galaxies and the 
general population of galaxies with SMBH in the universe. In 
the strictest sense the above results apply only to our sample 



of infrared selected AGN. However, if our sample is repre- 
sentative, which we reasonably expect it is, then it is also 
reasonable to extrapolate the above results to the general 
population of galaxies containing AGN. Furthermore, as the 
results of this analysis suggest that AGN do not appear to 
significantly distinguish between different galaxy morpholo- 
gies, at least at the level we can detect, it seems difficult to 
construct an argument to support the idea that the more 
general population of galaxies with a SMBH (i.e. active and 
inactive galactic nuclei) does. Moreover, if it is the case that 
the above statistical results can be extrapolated to the gen- 
eral population of galaxies containing a SMBH, they have 
an important corollary. Since it follows that if the current 
widely held opinion (and its supporting evidence) that most, 
if not all, massive ellipticals and galaxies with bulges con- 
tain a SMBH is valid, then essentially the same must apply 
to all the others, i.e. most galaxies do contain a SMBH. 

In comparing the redshift distribution Figure 4 with 
that of the galaxies in the SDSS DR 7 with structural pa- 
rameters, given in Figure 14 of Sll, we noted a sharper fall- 
off toward lower redshifts together with an associated shift 
in the peak of the distribution from about 0.1 to 0.2. This 
marked difference between the two distributions originates 
from the removal of apparently non-AGN galaxies from the 
sample at low redshifts by the WISE cross-match selection 
process. One possible reason for this apparent reduction of 
galaxies with AGN toward low redshifts could of course be 
the absence or paucity of SMBH in more recently formed 
galaxies. However, if this were the case, then as one would 
expect newly formed galaxies to be biased toward particu- 
lar morphology types, one should similarly expect to see a 



© 0000 RAS, MNRAS 000, 000-000 



8 Marleau et al. 




10^ 10"^ 108 1091010101110121013 10^ 10"^ 108 1091010101110121013 

Bulge Stellar Mass [M^] Disk Stellar Mass [Mq] 



Figure 9. Bulge {left) and disk stellar mass [right] versus blaclc liole mass obtained using the bolometric luminosity for our infrared 
sample of 5,718 galaxies not classified as bulgeless [black open circles). The bulge-dominated galaxies only {{B/T)r 5S 0.5) are highlighted 
as blue filled circles in the left figure and the disk-dominated galaxies only ({B/T)r < 0.5) are highlighted as blue filled circles in the 
right figure. As in Figure 8, the data points and bisector linear regression fit to the data (dashed line) are plotted for L^^l/ ^Edd = 0.1. 
Also shown are the fit for Li,oi/ LEdd = 1.0 {dotted line) and for Li,ai / LEdd = 0.01 {solid line). For comparison, we also plot the observed 
scatter of 0.3 dex of the empirical relation of Haring & Rix (2004) {grey shaded region). 



difference between the distribution of AGN in these mor- 
phology types with respect to the general population, which 
as we have seen is not obviously apparent. Hence, we are 
led to interpret the apparent depletion of AGN toward low 
redshifts as being more likely primarily due to the decrease 
in AGN activity toward low redshifts. Or, in other words, 
it would seem that the difference between the two redshift 
distributions should not be taken to imply that there exist 
fewer supermassive black holes in galaxies at lower redshifts, 
but merely that black holes become on average less active. 
This property, which is well known, can be seen for instance 
in the photometric redshift distribution of the SDSS quasars 
and their WISE counterparts, as depicted in Figure 2 of 
(Wu et al. 2012). Although here, by contrast, we stress that 
we are witnessing this effect in both the type 1 and type 2 
AGN galaxy populations. 



4.3 Stellar Mass Estimates 

Finally, we turn our attention to the relationship between 
the stellar and black hole masses for our infrared selected 
sample of galaxies with AGN. Following the description in 
Section 3.3 above, we obtained the stellar masses of 6,247 
galaxies in this sample from Mendel et al. (2012). In a sim- 
ilar fashion to our X-ray/radio sample, we used the Sersic 
profile total stellar mass estimates for our 530 bulgeless 
galaxies and the de Vaucouleurs -I- exponential profile total, 
bulge and disk stellar mass estimates for the remainder of 
our sample. The distribution of total stellar masses is shown 
in Figure 7 (left). 



4.4 Black Hole Mass Estimates from the 
Bolometric Luminosity 

Black hole masses were estimated using the bolometric lu- 
minosity of the galaxies. Bolometric luminosities (taken 
to be the 100 micron to 10 keV integrated luminosity 
Richards et al. 2006) are typically obtained using correc- 
tions to the mid-infrared bands where the AGN emission 
dominates. We used the 12 (W3) and 22 micron (W4) k- 
corrected flux densities from WISE to compute the bolo- 
metric luminosities of our WISE color-selected AGN and ap- 
plied the bolometric corrections Lboi — 8 x Li2iim (W3) and 
Lboi ~ 10 X L22Mm (W4) from Richards et al. (2006), which 
are not strongly dependent on AGN luminosity (see their 
Figure 12). Given the bolometric luminosity, making the 
assumption that accretion is at the Eddington limit yields 
a lower limit on the black hole mass. However, given that 
the growth rates of AGN can span several orders of mag- 
nitude, from super-Eddington accretion to W^'^LEdd, we 
computed black hole masses assuming the following three 
cases: L = LEdd, L = Q.lLsdd, and L = Q.QlLEdd- The 
histogram of BH masses (L = LEdd) is shown in Figure 7 
(right). The black hole masses range from 10^ Mq to 10* Mq 
with a peak near ~ 4 x 10®. These numbers are consistent 
with those measured in the sample of Simmons et al. (2011) 
(lO'^ Mq to 10*^ Mq) given our larger sample size. However, 
they are statistically lower than the masses measured from 
our X-ray/radio sample. 

Using the stellar masses of Mendel et al. (2012) and the 
black hole masses estimated from the bolometric luminosity 
of the galaxies, we present the correlation between black 
hole mass and total stellar mass for our sample of 6,247 
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galaxies, including the 530 classified as bulgeless in Figure 8 
(left). Here, we only show the results based on the 22 micron 
luminosities, although the results based on the 12 micron 
luminosities are very similar. 

The first thing to note is that the results of using our 
new method of selecting AGN based on their infrared col- 
ors clearly confirm the strong correlation between the total 
stellar mass of galaxies and the mass of their SMBHs. How- 
ever importantly, whereas previous studies were limited to 
primarily bulge-dominated systems, our study confirms this 
relationship to all morphological types, in particular, to bul- 
geless galaxies (Figure 8 left, green circles). Additionally, we 
emphasize that this study is by far the largest sample to 
date. 

In this figure, the total stellar mass estimates for the 
bulgeless galaxies are from the Sersic profile fitting (Figure 8 
left, green filled circles) whereas the rest are from the de 
Vaucouleurs bulge -I- exponential disk fitting (Figure 8 left, 
black open circles). For reference to previous work done for 
classical bulges, we show in the same figure the sample of 
galaxies with bulge Sersic index 3.9 < rib < 4.1 (Figure 8 
right, red filled circles). It may be noted that the location 
of these galaxies in the diagram is clearly indistinguishable 
from the bulgeless galaxies. 

A bisector linear regression fit to the data leads to the 
relation: 



log 



= a + j3 log 



IQii Mq 



(2) 



where a = 7.98 ± 0.20 and /3 = 1.14 ± 0.10 for 
Lboi/LEdd = 0.1. This relation is plotted in Figure 8 with a 
dashed line. The relations obtained using Lboi/LEdd ~ 1-0 
and Lfjoi/LEdd = 0.01 are also shown for comparison with a 
dotted and solid line, respectively. 

For comparison, we have also shown in the same fig- 
ure the empirical relation of Haring & Rix (2004) (where 
a = 8.20 ±0.10, /3 = 1.12 ±0.06, and scatter of 0.3 dex, grey 
shaded region). Although the slopes of the two fits are es- 
sentially the same, the agreement on the intercept depends 
on the value of Lboi/ LEdd- Assuming that previous studies 
such as that of Haring & Rix (2004) based on local dynam- 
ical mass estimates are correct, then our results appear to 
be most consistent with the black hole masses derived using 
Lboi/LEdd = 0.1. In fact, we expect this assumption to be 
the main source of uncertainty in our analysis, as the errors 
that arise from estimating the total stellar masses are on av- 
erage ~0.15 dex and those from estimating the bolometric 
luminosity are on average ~0.20 dex, and hence cannot ac- 
count for such a shift. Indeed, in this regard, Simmons et al. 
(2011) point out that accretion may take place at well below 
the Eddington limit and find that, at least for two sources in 
their sample, Lboi/LEdd = 0.05 and 0.08. Hence, if a growth 
rate of order Lboi/LEdd = 0.1 is typical of the sample, it 
implies that our black hole and total stellar mass relation 
is consistent with the empirical relation of Haring & Rix 
(2004). 

In Figure 9, we have plotted the results which essen- 
tially represent the galaxies in our sample with a significant 
bulge component, specifically those that do not satisfy the 
bulgeless selection criteria rig < 1.5 and {B/T)r < 0.1. In 
Figure 9 (left) we consider the bulge stellar mass only, while 
in Figure 9 (right) we consider the disk stellar mass only. In 
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Figure 10. Same as Figure 8, with galaxies divided in three 
redshift bins: z < 0.1 (blue), 0.1 ^ z < 0.2 (green), and z > 0.2 
(red). 



both figures scatter is significantly increased with respect to 
Figure 8, which uses the total stellar mass. This indicates 
that the true correlation is between the black hole mass and 
the total stellar mass. This is confirmed when we consider 
the bulge-dominated galaxies only {B/T)r ^ 0.5) in Figure 9 
(left, blue filled circles), and the disk-dominated galaxies 
only {B/T)r < 0.5) (right, blue filled circles), and note that 
these exhibit significantly tighter correlations. Therefore, we 
conclude that the previous assumption that the black hole 
mass was correlated with the bulge mass was only approxi- 
mately correct, since the true relation is in fact with the total 
stellar mass of the galaxy. However, the approximate bulge 
mass correlation therefore becomes increasingly accurate as 
the bulge mass dominates the total mass of the galaxy. 

In Figure 10, we also provide a plot of the data sorted 
into three redshift bins: z < 0.1 (blue open circles), 0.1 ^ 
2 < 0.2 (green open circles), and z ^ 0.2 (red open circles). 
Comparing data at low redshift, it is entirely consistent with 
what is found based on local dynamical mass estimates. In 
this figure, there appears to be a deficit of low masses at 
high redshift. Though we note that this could be purely due 
to the bias introduced by the limiting magnitudes of both 
the SDSS and WISE all-sky survey. 



5 SUMMARY 

By using the quantitative morphological analysis of 1.12 mil- 
lion galaxies in the SDSS DR7 (Sll), we have presented the 
first systematic survey of SMBHs in galaxies of all morpho- 
logical types, i.e. including bulgeless galaxies. Two different 
methods were used to identify the presence of black holes via 
their activity, the first being their X-ray/radio emission and 
the second their mid-infrared colors. These two studies are 
complementary in that while our first analysis of the X-ray 
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selected sample is biased towards selecting unobscured AGN 
and provides a well-tested and robust estimate of black hole 
mEiss but with a small statistic, our second analysis of the 
infrared selected sample includes both unobscured (type 1) 
and obscured (type 2) AGN and provides estimates on the 
black hole mass but with a large statistic. 

Our first analysis identified 688 unobscured active 
galaxies with known structural parameters. When the Sersic 
profile index Ug and bulge fraction {B/T)r of these galaxies 
were examined, we were able to identify 26 of these galax- 
ies hosting a SMBH as bulgeless. We were able to obtain 
estimates for both the stellar and black hole masses for 8 
of these 688 galaxies and investigate their relationship. The 
results of this first analysis confirm the correlation between 
black hole and total stellar mass for these 8 galaxies and 
includes one galaxy classified as bulgeless. 

Our second analysis identified f5,991 mostly obscured 
active galaxies with known structural parameters. When the 
Sersic profile index and bulge fraction of this sample were 
examined, we were able to identify f ,450 bulgeless galaxies 
containing AGN. Moreover, when analyzing the full range 
of structural parameters, we found that these active galax- 
ies followed a very similar morphological distribution to the 
general population of galaxies in the same redshift range. 
In particular, the fraction of bulgeless galaxies, ~ 9%, was 
found to be the same as in the general population. For this 
larger sample, we were also able to obtain estimates of the 
total stellar and black hole masses for 6,247 galaxies and in- 
vestigate their relationship. Importantly, whereas previous 
studies were limited to primarily bulge-dominated systems, 
our study confirms this relationship to all morphological 
types, in particular, to 530 bulgeless galaxies. Furthermore, 
we were able to investigate the same relationship using the 
independently estimated bulge and disk stellar masses. Our 
results indicated that the true correlation is between the 
black hole mass and the total stellar mass and hence we 
concluded that the previous assumption that the black hole 
mass is correlated with the bulge mass was only approxi- 
mately correct. If extrapolated, the combination of these re- 
sults appear to indicate that SMBHs are generic to galaxies 
regardless of the details of stellar galaxy evolution, dynam- 
ics and merger history of galaxies. As it is only possible to 
estimate dynamical masses locally, the future of studying 
SMBHs in galaxies at all redshifts lies in a method such as 
the one presented in this work. 
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